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Abstract 

I Adomian decomposition method is used for solving the seventh order 

' boundary value problems. The approximate solutions of the problems 

• I are calculated in the form of a rapid convergent series and not at grid 

'•pi ■ points. Two numerical examples have been considered to illustrate the 

I efficiency and implementation of the method. 
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The boundary value problems of ordinary differential equations play an important role in 
many fields. The theory of seventh order boundary value problems is seldom in the numer- 
ical analysis literature. These problems generally arise in modelling induction motors with 
two rotor circuits. The induction motor behavior is represented by a fifth order differential 
equation model. This model contains two stator state variables, two rotor state variables 
and one shaft speed. Normally, two more variables must be added to account for the effects 
of a second rotor circuit representing deep bars, a starting cage or rotor distributed pa- 
rameters. To avoid the computational burden of additional state variables when additional 
rotor circuits are required, model is often limited to the fifth order and rotor impedance 
is algebraically altered as function of rotor speed. This is done under the assumption that 
the frequency of rotor currents depends on rotor speed. This approach is efficient for the 
steady state response with sinusoidal voltage, but it does not hold up during the transient 
conditions, when rotor frequency is not a single value. So, the behavior of such models show 
up in the seventh order boundary value problems [lOj . 

Presently, the literature on the numerical solutions of seventh order boundary value prob- 
lems and associated eigen value problems is seldom. Siddiqi and Ghazala [III [131 [HI [13 
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[m [T71 [18] presented the solutions of fifth, sixth, eighth, tenth and twelfth order boundary 
value problems using polynomial and non-polynomial spline techniques. Siddiqi and Twiz- 
ell [121 [201 [21] presented the solution of eighth, tenth and twelfth order boundary value 
problems using eighth, tenth and twelfth degree splines respectively. 

The Adomian decomposition method was introduced and developed by George Adomian 
[H [21 [31 [3] . This method has been applied to a wide class of linear and nonlinear ordinary dif- 
ferential equations, partial differential equations, integral equations and integro-differential 
equations [HI [9l [221 [23] ■ Wazwaz [231 123] provided the solution of fifth order and sixth 
order boundary value problems by the modified decomposition method. Convergence of 
Adomian decomposition method was studied by Cherruault et al. [5l|6l[7]. In this paper, 
the Adomian decomposition method [11 [21 [3l [4] is applied on the seventh order boundary 
value problem. An algorithm for approximate solution by Adomian method is developed 
for such problems. This method provides solutions in terms of a rapidly convergent series. 



2 Adomian Decomposition Method 

Consider the operator form of differential equation 

Lu + Ru + Nu = f, (2.1) 

where L is the highest order derivative operator with the assumption that exists, R 
is a linear derivative operator of order less than L, Nu is the nonlinear term and / is the 
source term. Applying the inverse operator to both sides of Eq. (j 2. 1 p and using the 
given conditions yields the following 

u = g - L-^{Ru) - L-^{Nu), (2.2) 

where the function g represents the term obtained after integrating the source term / and 
using the given conditions . The solution u{x) can thus, be defined as 

oo 

"(^) = "^riix), ( 2. 3 ) 

n=0 

where the components uo,ui,U2, are determined by the following recursive relation 

uo{x) = g, 

Ufc+i(x) = -L-^R{uk) - L-^N{uk), k>0. ( 2. 4 ) 

The nonlinear term Nu is expressed in terms of an infinite series of the Adomian polynomials 
given by 

oo 

Nu = Y,An, ( 2. 5 ) 

n=0 

where An are Adomian polynomials [3j defined by 

, k = 0,1,2,... ( 2. 6 ) 

p=0 



An 



1 

n! dp" 
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3 Analysis of the Method 

Consider the following seventh order boundary value problem 

=(/>(x)n + ^/;(x) + /(2;,n), < x < b, (3.1) 

with boundary conditions 



itW(o) = ai, i = 0,1,2,3, 
u^'\b) = /3j, j = 0,1,2. 



(3. 2) 



where a^, i = 0,1,2,3 and Pj, j = 0,1,2 are finite real constants and the function / is 
continuous on [0,6]. 

The problem () 3. 1 p . in operator form, can be rewritten as 

Lu = (j){x)u + iIj{x) + f{x, u), ( 3. 3 ) 

where the derivative operator L is given by 

The inverse operator L^^ is therefore defined as 

rx px px px px px px 

= / / / / / / I {.) dx dx dx dx dx dx dx. ( 3. 5 ) 

Jo Jo Jo Jo Jo Jo Jo 

Operating with on Eq. () 3. 3 p gives 

u{x) = ao + aix + -^a2X^ + -^.a^x^ + -^a^x^ + -^a5X^ + -^oqx^ 

2\ 3! 4! 5! o! 

+L-\ct){x)u) + L-\i;{x)) + L-\f{x,u)), (3. 6) 

using the boundary conditions (| 3. 2 p at x = yields 

u(x) = ao+ aix + -^a2x'^ + -^asx^ + -^Ax"^ + -^Bx^ + ^Cx^ + 

2! 3! 4! 5! b! 

+L-\cl){x)u) + L-\^{x)) + L-\f{x, u)), ( 3. 7 ) 

where the constants 

^ = u(^)(0), B = u^^\0), C = n(6)(0), (3.8) 
will be determined using boundary conditions at j; = 6. 

The Adomian method determines the solution u{x) in terms of the following decomposition 
series 

oo 

"(^) = XI ^n{x). ( 3. 9 ) 

n=0 
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Substituting Eq. (| 3. 9 p into Eq. (| 3. 7 p leads to 

°° 11111 
^u„(x) = ao+ aix + —a2X^ + —a^x^ + ^Ax^ + —Bx^ + -^Cx^ + 

n=0 

oo 

+L-\(t>{x) J2 M^)) + L~\^{x)) + L~\f{x, u)). ( 3. 10 ) 

n=0 

To determine the components Un{x), n >0, the fohowing recurrence relation will be used 

uo(x) = ao + aix + -^a2x'^ + -^asx^ + -^Ax"^ + -^Bx^ + ^Cx^ + L^^(ip(x)), 

2! 3! 4! 5! b! 

oo oo 

Uk+lix) = L-\(t>{x)Y,Uk{x)) + L~\Y.^k), A:>0, (3.11) 

fc=0 k=0 

where the nonlinear function f{x,u) can be decomposed into an infinite series of Adomian 
polynomials given by 

oo 

f{x,u) = Y,Ak, (3.12) 

A:=0 

where Ak are Adomian polynomials [3] defined by 

1 



oo 



, n = 0,1,2,... ( 3. 13 

p=0 



Keeping in view Eq. (j 3. 11 p - Eq. (j 3. 13 p . the components no,ni,U2,... are calculable. 
By the Adomian method the solution can be constructed as 

u = lim ifn, 

n— >oo 

where the n-term approximant is defined by 

n-l 
1=0 

Applying the boundary conditions at x = 6 to the approximant ipn- The resulting algebraic 
system in A, B and C can be solved to produce approximations to constants A, B and C. 
Finally, the approximated solution of the seventh order boundary value problem follows. 
To implement the method, four numerical examples are considered in the following section. 

4 Numerical Examples 

Example 4.1 Consider the linear seventh order boundary value problem 

= xu{x) + e''{x^ -2x-6), < x < 1, (4.1) 
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(4. 2) 



subject to the boundary conditions 

u{0) = 1, n(l) = 0, 

^,{i)(0) = 0, u«(l) = -e, 

u(2)(0) = -1, u(2)(l) = -2e, 

u(3)(0) = -2. 

The exact solution of the problem (4.1) is 

u{x) = (1 - 2;)e^. 

The problem (| 4. 1 p . in operator form, can be rewritten as 

Lu = xu{x) + e''{x'^ -2x -6), (4.3) 

Operating with L^^ on Eq. (j 4. 3 P and using the boundary conditions (j 4. 2 p at x = 
gives 

u{x) = -63-6ix-^x'^ - 4x^ + {-- + — )x^ + {-— + — )Bx^ + — {2 + C)x^ 
^ ^ 2! ^ 2 24^ ^ 30 120^ 360^ ^ 

+e^(-8 + xf + L^^{xu{x)), ( 4. 4 ) 

where the constants 

A = n(^) (0), B = (0), C = 'u(^) (0), (4.5) 

are to be determined. Substituting the decomposition series (j 2. 3 p for u{x) in Eq. (j 4. 4 p 
yields 

Yunix) = -63-64x-^x^-4x^ + {-- + — )x* + {-— + — )Bx^ + — (2 + C)x'^ 

2! ^2 24^ ^ 30 120^ 360^ ' 

oo 

+e^(-8 + xf + L~\x J2 Mx))- ( 4. 6 ) 

n=0 

Using the recurrence algorithm (| 3. 11 P the following relations are obtained 

Qc: 14 1 R 1 

uo(x) = -63 - 64x - - 4x3 + ( — + —)x^ + ( + )Bx^ + (2 + C7)x^ 

uw 2! 2 24^ ^ 30 120^ 360^ ^ 

+e"(-8 + x)2, 

Wfc+i(x) = L-^(xufc(x)), A; > 0. ( 4. 7 ) 

From Eq. (j 4. 7 p . 

no(x) = -63-64x-^x2-4x3 + (-i + ^)x4 + (-l + A)i53;5^_L(2 + C)x6 

+e"(-8 + x)2, 
ui{x) = L"^(xno(x)), 

= 1848 + 1392x + 505x2 ^ ^^g^3 _^ 4 _^ 5 _^ 



2 15 120 640 3780 34560 
xii (-12 + ^)xi2 (_4 + 5)2;i3 (2 + C7)xi4 

+ nrronnoon + \ no^oo^onn + /o.^.n.-i^nn + e''(-1848 + 456x 



415800 95800320 1037836800 12454041600 
-37x2 + x3), 

(4. 8) 
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The series solution of u{x) is found to be as an approximation with four components 

u{x) = uo{x) + ui{x) + U2{x) + U3{x) . (4.9) 

The unknown constants A, B and C can be obtained by imposing the boundary conditions 
at X = 1 on the four-term approximant (^4 defined by Eq. () 4. 9 p , as 

A = -3.0000001004600083, B = -3.9999991664171546, C = -5.000002140841064. 

Finahy, the series solution can be written as 



U[X) 



1 



23 7 
^ ^ 0.125x^ - 0.0333333x^ -0.00694445x^ ^ 



X 



X 



X 



2 3 
10 



403200 3991680 

14\ 



+0{x 



840 5760 45360 
(2.29644 X 10"^)x^2 - (1.60591 x 10"^°)x^^ 

( 4. 10 



The comparison of the exact solution with the series solution of the problem (4.1) is given in 
Table 1, which shows that the method is quite efficient. Figure 1 also endorses the efficiency 
of the method. 

Example 4.2 Consider the following seventh order nonlinear boundary value problem 



«(^)(x) 

subject to the boundary conditions 



e^u (x), < X < 1, 



4. 11 ) 



n(0) = 


1, 


n(i; 


uW(0) = 


-1, 


n«(l 


n(2)(0) = 


1, 


n(2)(r 


n(3)(0) = 


1 

1 — ^ 





The exact solution of the problem (4.2) is 

u{x) = e~^. 

The problem () 4. 11 p . in operator form, can be rewritten as 

Lu = e^n^(x). 



( 4. 12 ) 



4. 13 ) 



Operating with L ^ on Eq. () 4. 13 p and using the boundary conditions () 4. 12 p at x = 
gives 

u{x) = 1 - X + ^x2 - i-x^ + ^Ax"^ + ^Bx^ + ^Cx^ + L-i(-e^n2(x)), ( 4. 14 ) 
2! 3! 4! 5! 6! 

where the constants 



4. 15 ) 
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Figure 1 : Comparison of the approximate 
solution with the exact solution for the 
problem (4.1). Dotted line: approximate 
solution, solid line: the exact solution. 




Figure 2: Comparison of the approximate 
solution with the exact solution for the 
problem (4.2). Dotted line: approximate 
solution, solid line: the exact solution. 



are to be determined. Substituting the decomposition series (j 2. 3 p for u{x) and the series 
of polynomials (I 2. 5 p for u^(x) in Eq. (14. 14 p yields 



n=0 



+ L-i(-e^E^")' (4-16) 



n=0 

,2/ 



where An are the Adomian polynomials for the nonlinear term Nu = u (x), thus given by 

^0 = N{uo) = ul{x), 

Ai = uo{x)N (uq) = 2uo{x)ui{x), 



2/ \ 

A2 = U2{x)n'{uq) + ^^iV"(no) = 2uo(x)u2(x) + u\{x), 



( 4. 17 ) 



Using the recurrence algorithm (| 3. 11 p yields 

uq{x) = 1 - X + - + Iax^ + isx^ + 
Uk+i{x) = L-\-e'Ak), k>0, (4.18 
Prom (| 4. 17 p and (| 4. 18 p . 

no(x) = 1 - X + ^x"^ - ^x^ + ^Ax'^ + ^Bx^ + ^Cx*^, 
^ ^ 2! 3! 4! 5! 6! 

ni(x) = L-i(-e^Ao), 



x'^ X* x*^ x^° , 1 A , 
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5040 40320 362880 3628800 '39916800 19958400' 



1 + f )a;i2 + ( J - -^^y^ + 0{x'^) 



^479001600 239500800' '6227020800 3113510400' 
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Figure 3: Absolute Error for the problem Figure 4: Absolute Error for the problem 
(4.3). (4.4). 



The series solution of u{x) is found to be as an approximation with two components 

u{x) = uo{x) + ui{x). (4. 19) 

The unknown constant A, B, and C can be obtained by imposing the boundary conditions 
at X = 1 on the two-term approximant ip2 defined by Eq. (14. 19 P , as 

A = 1.0000000197456873, B = -1.000000228700079, C = 1.0000008112153262. 

Finally, the series solution can be written as 

u{x) = 1-X + 0.5x2 _ o.i66667x^ + 0.0416667x^ - 0.00833334x^ + 0.00138889x^ 

-0.000198413x^ + 0.0000248016x^ - (2.75573 x lO"^)^^ + (2.75573 x 10"^)x^° 

-(2.50521 X 10~®)x^^ + (2.08768 x lO"^)^^^ - (1.60591 x 10"^°)x^^ 

+0{x^^). ( 4. 20 ) 

The comparison of the exact solution with the solution of the problem (4.2) is given in Table 
2, which shows that the method is quite efficient. Figure 2 also endorses the efficiency of 
the method. 



Example 4.3 Consider the following seventh order linear boundary value problem 

< X < 1, 



uC^)(x) = -u(x) -e^(35 + 12x + 2x2) 

u{0) = 0, n(l)=0, 

uW(0) = 1, u«(l) = -e, 

(2) (0) = 0, n(2)(l) = -4e, 

(3) (0) = -3. 



u 



u 



( 4. 21 



The exact solution of the Example 3.1 is u{x) = x(l — x)e^ [11] . 

Following the procedure of the Example 4.1, this problem is solved. It is observed that the 
errors in absolute values are better than those of Siddiqi and Iftikhar [11] as shown in Table 
3. In Figure 3 absolute errors are plotted. 



Example 4.4 The following seventh order nonlinear boundary value problem is considered 
= u(x)'u'(x) + e-2^(2 + e^(x-8) -3x + x2), < a; < 1, ^ 



u{0) = 1, u{l) = 0, 

„(i)(0) = -2, uW(l) = -l/e, 

u(2)(0) = 3, u^'^\l) = 2/e, 

u(3)(0) = -4. 



( 4. 22 



The exact solution of the Example 3.2 is u{x) = (1 — x)e~^. 

Following the procedure of the previous problem 4.2, this problem is solved. The comparison 
of the exact solution with the solution of the problem 4.2 is given in Table 4, which shows 
that the method is quite efficient. In Figure 4 absolute errors are plotted. 
Conclusion In this paper, the Adomian decomposition method has been applied to obtain 
the numerical solutions of linear and nonlinear seventh order boundary value problems. The 
numerical results show that the method is quite efficient for solving high order boundary 
value problems arising in various fields of engineering and science. 



Table 1: Comparison of numerical results for t 


le problem 4.1 


X 


Exact solution 


Approximate series solution 


Absolute Error 


0.0 


1.0000 


1.0000 


0.0000 


0.1 


0.0994 


0.0994 


4.3972E-10 


0.2 


0.9771 


0.9771 


4.9251E-10 


0.3 


0.9449 


0.9449 


7.4067E-10 


0.4 


0.8950 


0.8950 


6.6537E-10 


0.5 


0.8243 


0.8243 


3.0059E-11 


0.6 


0.7288 


0.7288 


4.3591E-10 


0.7 


0.6041 


0.6041 


3.6735E-10 


0.8 


0.4451 


0.4451 


7.2753E-10 


0.9 


0.2459 


0.2459 


7.0036E-10 


1.0 


0.0000 


2.2191E-10 


2.2191E-10 



Table 2: Comparison of numerical results for problem 4.2 



X 


Exact solution 


Approximate series solution 


Absolute Error 


0.0 


1.0000 


1.0000 


0.0000 


0.1 


0.9048 


0.9048 


1.5676E-9 


0.2 


0.8187 


0.8187 


1.6418E-9 


0.3 


0.7408 


0.7408 


4.9680E-9 


0.4 


0.6703 


0.6703 


1.5514E-9 


0.5 


0.G065 


().()()()3 


1.5274E-9 


0.6 


0.5488 


0.5488 


2.4958E-9 


0.7 


0.4965 


0.4965 


1.3993E-8 


0.8 


0.4493 


0.4493 


2.5593E-9 


0.9 


0.4065 


0.4065 


5.4089E-9 


1.0 


0.3678 


0.3678 


1.1034E-9 
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Table 3: Comparison of numerical results for the problem 4.3 



X 






iT-UoUlLlljC J_Jii(Ji 
LJJ. ^jO'^IJ.Ij J.1J.O Li J.1WV4. 


Siddini and Tftikhar 111 


0.0 


0000 


0000 


0000 


0000 


0.1 


0.9946 


0.9946 


1.23082E-13 


8.55607E-13 


0.2 


0.1954 


0.1954 


3.7792E-13 


9.94041E-12 


0.3 


0.2835 


0.2835 


2.37421E-13 


3.52244E-11 


0.4 


0.3580 


0.3580 


3.62099E-13 


7.3224E-10 


0.5 


0.4122 


0.4122 


9.39249E-14 


1.08769E-10 


0.6 


0.4373 


0.4373 


4.82947E-13 


1.29035E-10 


0.7 


0.4229 


0.4229 


1.09135E-13 


1.51466E-10 


0.8 


0.3561 


0.3561 


1.64868E-14 


2.717974E-10 


0.9 


0.2214 


0.2214 


7.25975E-13 


7.48179E-10 


1.0 


0.0000 


-4.54747E-13 


4.54747E-13 


2.1729E-09 



Table 4: Comparison of numerical results for the problem 4.4 



X 


Exact solution 


Approximate series solution 


Absolute Error 


0.0 


1.0000 


1.0000 


1.67932E-12 


0.1 


0.814354 


1.1051 


2.96696E-12 


0.2 


0.654985 


0.654985 


1.26055E-12 


0.3 


0.518573 


0.518573 


2.10898E-12 


0.4 


0.402192 


0.402192 


6.68926E-12 


0.5 


0.303265 


0.303265 


7.21923E-12 


0.6 


0.219525 


0.219525 


9.75339E-12 


0.7 


0.148976 


0.148976 


2.19552E-12 


0.8 


0.0898658 


0.0898658 


4.24917E-12 


0.9 


0.040657 


0.040657 


2.27311E-13 


1.0 


0.0000 


4.42298E-12 


4.42298E-12 
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